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. We are interested in understanding the dynamics of dissipative partial differential equa- 

tions on unbounded spatial domains. We consider systems for which the energy density e > 
satisfies an evolution law of the form dte — div^f ~ d, where —/is the energy flux and d> 
the energy dissipation rate. We also suppose that |/p < b{e)d for some nonnegative function 
6. Under these assumptions we establish simple and universal bounds on the time-integrated 
fyf) ' energy flux, which in turn allow us to estimate the amount of energy that is dissipated in a 

. given domain over a long interval of time. In low space dimensions N < 2, we deduce that 

' any relatively compact trajectory converges on average to the set of equilibria, in a sense 

that we quantify precisely. As an application, we consider the incompressible Navier-Stokes 
CN| ■ equation in the infinite cylinder M x T, and for solutions that are merely bounded we prove 

that the vorticity converges uniformly to zero on large subdomains, if we disregard a small 
subset of the time interval. 



1 Introduction 

Many time-dependent partial differential equations arising in Mathematical Physics are dissipa- 
tive in the sense that there exists a nonnegative energy density e{x,t), depending on the space 
variable x G and the time t, which is locally dissipated under the evolution defined by the 
system. By this we mean that e{x,t) satisfies an equation of the form 

dte{x,t) = dW^f{x,t) - d{x,t) , (1) 

where —f{x,t) G denotes the energy flux in the system and d{x,t) > is the energy 
dissipation rate. Equivalently, integrating ([1]) with respect to both variables x, t and applying 
the divergence theorem, we obtain the energy balance equation 

/ e{x,T2)dx- / e{x,Ti)dx = / / f{x,t)-udadt- / / d{x,t)dxdt, (2) 
Jn Jn Jti Jdn Jti Jn 

which holds for all Ti < T2 and all admissible domain Q C M^. Here u denotes the outward 
pointing unit normal on dO,, and da is the elementary surface area. 
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As a typical example, consider the reaction-diffusion equation 

dtu{x,t) = A-,u{x,t) -V'{u{x,t)) , xeR^, t>0, (3) 

where and y : R — )• M+ is a smooth potential. This equation appears for 

instance in the theory of phase transitions [3] and in population genetics |5j. In the particular 
case where V{u) = |(1 — u^)"^, Eq. ([3]) is often referred to as the Allen-Cahn equation or the 
real Ginzburg-Landau equation. If u{x,t) is any smooth solution of we define the energy 
density, the (backward) energy flux, and the energy dissipation rate by the formulas 

e = ^|Vn|2 + V{u) , / = utVu , d = , (4) 

where ut = dfU. It is then straightforward to verify that ([T]) holds, which means that energy is 
locally dissipated under the evolution defined by ([3]). Since V is nonnegative, we also deduce 
from (j4]) that |/p < 2ed. We shall list in Section [2] several other examples of classical PDE's 
which define dissipative dynamical systems in the same sense. In most of these examples, the 
energy flux happens to satisfy an inequality of the form 

I/I' < Ced , (5) 

for some positive constant C. 

If a dissipative PDE such as ^ is considered in a bounded domain C M^, with boundary 
conditions ensuring that / • < on dil., then ([2]) shows that the total energy 

E{t) = / e{x,t)dx 
Jn 

is a Lyapunov function of the system, namely E{t) is a decreasing function of time for all 
solutions of ([3]) which are not equilibria. Under natural coercivity assumptions on the potential 
V, this gradient structure implies that all finite-energy solutions of (l3|) in a bounded domain 
converge to the set of equilibria as t — )• oo [I9l [20] . 

The situation is rather different if we work in an unbounded domain such as = M^. In 
that case, Eq. ([3]) may have travelling wave solutions of the form u{x,t) = v{x — ct) for some 
nonzero c G [5], and such solutions do not converge uniformly to equilibria as t — )■ oo. One 
may object that travelling waves do converge to equilibria uniformly on compact sets, but it is 
possible to construct more complicated solutions for which convergence to equilibria does not 
hold even in that weaker sense, see Example 15.71 below. Moreover, if > 2, one can exhibit 
scalar reaction-diffusion equations of the form dtu = Au + F{x, u) which have nontrivial time- 
periodic solutions |15j . This is in sharp contrast with what happens for gradient systems, but 
one should keep in mind that all counter-examples above involve infinite-energy solutions. 

When the total energy cannot be used as a Lyapunov function, a natural idea is to exploit 
the energy balance equation ([2|) or its differential version ([1]) to obtain relevant information 
on the dynamics. In the context of extended dissipative systems, this approach was initiated 
in a previous paper by the authors [15], the main conclusions of which can be summarized as 
follows. If A < 2, the reaction-diffusion equation ([3]) on cannot have any nontrivial solution 
such that e{x,T2) > e{x,Ti) for some T2 > Ti and all x G M^; in particular, nontrivial time- 
periodic solutions are excluded. Furthermore, all bounded solutions converge on average (in 
time), uniformly on compact sets (in space), toward the set of equilibria as t — )• +00. In other 
words, due to the local energy dissipation law ([1]), Eq. ([3]) retains some dynamical properties of 
usual gradient systems, provided A < 2. In contrast, if A > 2, highly non-gradient behaviors 
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such as nontrivial periodic orbits can occur. The conclusions of [T5] also apply to a damped 
hyperbolic equation which reduces to ([3]) in the limit of strong damping. 

The aim of the present paper is to strengthen and generalize the results of |15j . Instead of 
considering a particular equation, we work in an abstract setting, assuming only the local energy 
dissipation law ([T|) and an estimate of the form ([5|) for the energy flux. As a consequence, our 
results apply to a much larger class of systems, some of which are listed in Section [21 Another 
substantial progress with respect to [15j is a new estimate on the time-integrated energy flux 
through a closed hypersurface, which we derive in Section [3l This bound allows us to obtain 
quantitative versions of the main results in [15]. For instance, in Section U] we give an explicit 
estimate of the energy dissipated in a given domain over a long time interval, and in Section [5] we 
measure the fraction of time spent by any relatively compact trajectory outside a neighborhood 
of the set of equilibria. As can be expected from [15j, our results depend strongly on the space 
dimension A^, and some of them even fail if A^ > 2. As a final application, we consider in 
Section [6] the incompressible Navier-Stokes equation in the two-dimensional cylinder M x T, and 
for solutions that are merely bounded we prove some convergence results for the vorticity which 
are apparently new in this context. 

2 Extended Dissipative Systems 

To treat in a unified way various dissipative PDE's on unbounded domains, we introduce in this 
section the notion of an extended dissipative system, which will be studied in the rest of this 
paper. We also list a few classical examples which fit into our abstract framework. 

Let A be a metrizable topological space. We say that a family {^(t))t>o of continuous maps 
in A is a continuous semiflow on X if 

i) <I>(0) = 1 (the identity map); 

ii) + t2) = $(ti) o$(t2) for all ti,t2 > 0; 

iii) For any T > 0, the map {t,u) i— )• ^{t)u is continuous from [0,T] x A to A. 

In particular, if uq G A, the trajectory u : M+ — t- A defined by u{t) = ^{t)uQ for all 
t G M+ = [0,oo) is continuous, and u{t) depends continuously on the initial data uq, uniformly 
in time on compact intervals. As an example, if V{u) = |(1 — u^)^, the reaction-diffusion 
equation ([3D defines a continuous evolution semifiow on the space A = C'^(M^) n L°°(]R^), if A 
is equipped with the topology of uniform convergence on compact sets of M^. The systems we 
are interested in are those for which one can define an energy density e, an energy fiux /, and 
an energy dissipation rate d with the same properties as in the example This leads to the 
following definition : 

Definition 2.1. Let A € N*. We say that a continuous semiflow {^{t))t>o on a metrizable 
space X is an extended dissipative system on M.^ if one can associate to each u £ X a triple 
{ej,d) with e,de C°{R^,R+) and f G C°(M^,M^) such that: 

Al : The functions e,f,d depend continuously on u £ X, uniformly on compact sets o/M^; 
A-2: l/p < b{e)d for some nondecreasing function b : M+ — )• M+; 
A3: d = only if^{t)u = u for all t > 0; 

A4 ■' Under the evolution defined by the semiflow {^{t))t>o, the time- dependent quantities e, f, d 
satisfy the energy balance dte = div/ — d in the sense of distributions on M.^ x M.-^.. 

Remarks 2.2. 

1. More generally, one can define extended dissipative systems on any (unbounded) domain 



3 



$7 C by substituting Q for everywhere in Definition 12.11 For simplicity, we restrict 
ourselves to the case of the whole space in this paper, except in Section [6] where we shall 
consider a system defined in a cylindrical domain. 

2. We emphasize that, in Definition 12.11 both the energy density e and the energy dissipation 
rate d are supposed to be nonnegative. The first condition ensures that the energy density is 
bounded from below, and the positivity of d together with the energy balance ([T]) imply that 
energy is locally dissipated (and never created) in the system. 

3. In assumption A2, it is understood that the function b : — )• M+ is independent of 
u & X. In many examples one can take 6(e) = Ce for some positive constant C, as in ([5]), 
but the generalization proposed here is necessary if one considers systems such as the nonlinear 
diffusion equation (see below) or the two-dimensional vorticity equation (see Section [6]) . 

4. Assumption A3 means that all trajectories, except equilibria, dissipate some energy. Note 
that, if d = 0, then / = by A2, hence ([T]) already implies that the energy density is time- 
independent. We emphasize that an extended dissipative system may have equilibria for which 
d is not identically zero; these may be called nonequilibrium steady states, in the terminology 
of Statistical Mechanics. On the other hand, if one considers systems with a continuous group 
of symmetries, it may be useful to relax assumption A3 so as to allow for a vanishing energy 
dissipation on relative equilibria of the system; these are equilibria up to the action of the 
symmetry group, see the example of the complex Ginzburg-Landau equation below. 

5. Our final comment on Definition 12.11 concerns the regularity of e, /, and d. To avoid techni- 
calities, we have supposed that, for each u £ X, the densities e, /, d are continuous functions on 
M^. Moreover, if u{t) varies continuously in time, the associated quantities e{t) , f (t) , d{t) are 
jointly continuous in space and time, in view of assumption Al. In the particular case where 
u{t) = ^{t)uQ for some uq G X, assumption A4 implies that the integrated energy balance 
equation ([2]) holds for all T2 > Ti > and any smooth bounded domain Q C M^, and since $ 
is a continuous semiflow we even know that all four terms in ([2]) depend continuously on the 
initial data uq £ X. These comfortable assumptions are not unrealistic, and can be verified in 
all systems listed below if we choose functions spaces of sufficiently high regularity. However, 
especially in nonparabolic PDE's, it is often more convenient to use larger function spaces, in 
which (for instance) the energy density is locally integrable but not continuous. In that case, 
intead of Al and A4, it is sufficient to require that the energy balance equation ([2]) be satisfied 
for all r2 > Ti > and any smooth bounded domain C M^, and that the various quantities 
in ([2]) depend continuously on the initial data. 

In complement to Definition 12.11 we specify two properties of extended dissipative systems 
which will play an important role in Sections [3] to [H 

Definition 2.3. We say that an extended dissipative system in the sense of De finition \2. i\ is 

• bounded, if there exists E > such that e{x) < E for all x G and all u G X; 

• compact, if the space X is compact. 

Remark 2.4. In many extended dissipative systems, boundedness can be achieved simply by 
restricting the space X to a subset which is positively invariant under the evolution defined by 
the semiflow, and on which the energy density is uniformly bounded by some positive constant 
E. In that case, we can replace 6(e) by a positive constant f3 > b{E) in assumption A2. We 
thus obtain the relation d > which quantifles how much energy is dissipated in the 

system when the energy flux is not identically zero. On the other hand, compactness can often 
be achieved by endowing X with a sufficiently weak topology. Note however that, with the 
definitions above, a compact extended dissipative system is not necessarily bounded. 
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To substantiate Definitions 12.11 and 12.31 we now give a few concrete examples. 
Examples 2.5. 

1. A reaction-diffusion equation 

We consider again the reaction-diffusion equation ([3|), and specify in wliich function spaces it 
defines an extended dissipative system in tlie sense of Definition 12.11 There are of course many 
possibilities, and we just mention here two reasonable ones. Since we want global solutions of ([3]), 
it is natural to assume that the potential y : M — t- is coercive in some sense. For instance, we 
can suppose that there exists a constant m > such that uV'{u) > m whenever \u\ is sufficiently 
large. Then it is known that the Cauchy problem for ([3]) is globally well-posed in C^^(]R^), the 
Banach space of all functions u : — )• R that are bounded and uniformly continuous together 
with their derivatives up to order k £ N. This means that ([3]) defines a continuous semiflow 
{^{t))t>Q on X = C^^(R^). Moreover, if A; > 2, the quantities e,f,d introduced in ([4]) belong 
to C'^(R^) and depend continuously on it G X, uniformly on compact sets of R^. Together 
with ([I]) and dSI), this shows that, if /c > 2, the semiflow of ^ on X = C^^{R^) is an extended 
dissipative system. 

This system is not bounded in the sense of Definition 12.31 but it becomes so if we choose for 
X any bounded subset B C C^y(R^) that is positively invariant under the semiflow. In addition, 
if > 2 and if we consider i3 as a subset of C''(R^), equipped with the topology of uniform 
convergence on compact sets, then the closure of B in C'^(M^) is compact and the restriction 
of the semiflow {^{t))t>o to B defines a compact extended dissipative system. The idea of 
introducing a localized topology to restore compactness plays an important role in the study 
of dissipative PDE's on unbounded domains, in particular when constructing global attractors 

Instead of C^^(R''^), another possible choice is the uniformly local Sobolev space H^^iMJ^), 
on which ^ also defines a continuous semifiow if s > N/2, see [6l [15]. If moreover s > 2 + N/2, 
the densities ([4]) are continuous and we again obtain an extended dissipative system. As above, 
if we restrict our analysis to a bounded invariant subset B C //^^(R^) and if we take the closure 
of B with respect to the topology of Lf^^{MJ^), the restriction of the semifiow to B defines an 
extended dissipative system that is bounded and compact in the sense of Definition 12. 3i 

2. A strongly damped wave equation |24] 

Given a > and a smooth potential F : R — t- we consider the equation 

utt + ut - aAut = I^u - V {u) , xGR^, t>0, (6) 

where u : x R+ — )■ R. As usual, this second-order equation can be written as a first-order 
system for the pair (u,ut). For simplicity, we assume that the potential V is quadratic near 
infinity, namely V"{u) = m > for all sufficiently large u G R. Then the initial value problem 
for Eq. Q is globally well-posed in the uniformly local space X = H^yCR^) x i?^j"-'^(M^) if 
s > N/2. For any pair {u,ut) G X, we introduce the densities 

e = ^uf + ^|Vn|2 + V{u) , / = ut{Vu + aVut) , d = li? + a\Vut\'^ , 

which are well-defined and continuous provided s > 2 + A^/2, or s > 1 + A^/2 if a = 0. Then ([6]) 
implies that the local energy dissipation law ([1]) is satisfied, and a direct calculation shows that 
([U holds with C = 2max(l, a). Finally d = implies ut = 0. Thus the semifiow of the strongly 
damped wave equation ([6]) in X is an extended dissipative system in the sense of Definition 12. 11 
In the particular case where q = 0, the local dissipation of energy for Eq. ([6]) was already studied 
in [E]. 
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3. The complex Ginzburg-Landau equation [H [9l [22] 

Our next example originates from the complex Ginzburg-Landau equation 

ut = {l + ia)Au + u- {l + il3)\ufu , x G , t > , (7) 

where u : x IR+ — C and a, /3 are real parameters. To have a gradient structure, we assume 
that /3 = a, and we introduce the auxiliary function v{x,t) = w(x,i)e*"*, which satisfies the 
equation 

vt = {l + ia)(^Av + v -\v\'^v'j , x € , t > . (8) 

The Cauchy problem for ([8]) is globally well-posed in the function X = C^y(M^,C) for A; > 
or X = i?^i(R^, C) for s > N/2. If in addition A;>2ors>2 + N/2, then for any f G X the 
densities ^ 

e = h^v\' + ]{l-\v\^f , f = Re{vtVv), d = , 
2 4 1 + 

are well-defined and continuous. A direct calculation also shows that ([1]) holds, as well as ([5]) 
with C = 2{l + a^). Thus the semiflow of ([5]) in X is an extended dissipative system in the sense 
of Definition 12.11 We also remark that d = if and only if u{x,t) = ?;(x)e~*"* for some v £ X, 
which means that u{-,t) is a relative equilibrium of ([7]) : u{-,t) moves without dissipation along 
an orbit of the symmetry group U{1). Thus the semiflow of ([7]) in X is an extended dissipative 
system only if assumption A3 is relaxed as suggested in Remark 12.21 4. 

4. The Landau- Lifshitz- Gilbert equation |18j 

We now consider a vector- valued PDE appearing in micromagnetism. Given a G M, the Landau- 
Lifshitz equation reads 

ut = -uA{uA An) + an A An , x , t>0 , (9) 

where u : x — )• = {v G M^llt;! = 1}. Here A denotes the usual cross product in M^. 
In particular —u A (n A An) = An — n(n • An) = An + |Vnpn is the orthogonal projection of 
An onto the plane orthogonal to the direction n G S^, and n A An is the same vector rotated by 
7r/2 in the orthogonal plane. The initial value problem for ([9]) is locally well-posed in the space 
X = C^^(M^) for /c > or X = H^^{R^) for s > N/2, but in general finite-time singularities are 
expected to occur, unlike in the previous examples. To obtain a continuous semiflow {^{t))t>o, 
it is therefore necessary to restrict our space X to a family of global trajectories. Now, if A; > 2 
or s > 2 + X/2, the densities 

1 ^ 
e = 21^""!^ ' / = ""t^^ = '^(.dtUk)Vuk , d = |nAAn|^ , 

fe=i 

are well-defined and continuous for any n G X, and it is again straightforward to verify that ([I]) 
and ([5]) hold with C = 2(1 + a"^). Moreover d = implies n A An = 0, hence nj = 0. Thus 
Eq. ([9]) also defines an extended dissipative system in the sense of Definition I2.H provided we 
restrict the space X to a suitable family of global solutions. 

5. A nonlinear diffusion equation 

To motivate assumption A2 in Definition [2Tj we also give an example where the relation between 
the energy flux and the energy disipation is more complex than in ([5]). Given a smooth function 
a : M — )• (0, oo), we consider the nonlinear diffusion equation 

Ut = div(a(n)Vn) , x G , t>0 , (10) 



6 



which is globahy weh-posed in the space X = C^u(M^) for > or X = H^\{M.^) for s > N/2. 
IfA;>lors>l + N/2, we denote for all u G X : 

e = -u^ , / = ua{u)Vu, d = a(M)|VM|^ . 

Then ([T]) holds, and since a{u) > it is clear that d = implies Uf = 0. Moreover |/p < 2a{u)ed. 
Thus, if we define 

5(e) = 2esup{a(ii) <2e} , e > , 

then e i— )• 6(e) is increasing and |/p < 5(e) d by construction. Thus (jlOp defines an extended 
dissipative system in X in the sense of Definition I2.1i 

6. T/ie two-dimensional vorticity equation 

As a final example, we consider the vorticity equation associated to the two-dimensional incom- 
pressible Navier-Stokes system. In this model, the velocity of the fiuid, which is denoted by 
u{x,t) G M^, satisfies the incompressibility condition diui + d2U2 = 0, and the corresponding 
vorticity field oj = diU2 — d2Ui evolves according to the advection-diffusion equation 

Sfo; + n • VtJ = Aw , x G , t>0. (11) 

If we define the enstrophy density e, the enstrophy fiux /, and the enstrophy dissipation rate d 
by the formulas 

e = -uj^ , f = ujVuj — 2^^"^ ; d = |Va;p , (12) 

it is easy to verify that ([1]) is satisfied. Moreover, d = clearly implies that dfOJ = 0. However, it 
is not possible to obtain here an inequality of the form ([5]), nor of the generalized form adopted 
in assumption A2 of Definition 12.11 The main difficulty comes from the term ^uuj"^ in the 
enstrophy fiux. Since the velocity u is not a local function of to, the bound ([5]) cannot hold 
pointwise unless the contribution of u is absorbed into the constant C in the right-hand side. 
This requires a uniform bound on the -L°° norm of the velocity field, which is not known to hold 
for solutions of the two-dimensional Navier-Stokes equations that are only bounded in space, see 
\n\ [26l [27] . In addition, the term |ua;^ does not contain any derivative of oj, hence does not 
necessarily vanish when d = 0. This means that enstrophy can (a priori) be transported without 
any dissipation, whereas it is essential in our approach that the energy dissipation be bounded 
from below in terms of the energy flux. Surprisingly enough, we shall show in Section[6]that these 
difficulties essentially disappear if we consider the vorticity equation (|lip in the infinite cylinder 
R X T instead of the whole plane M?. Thus, if we assume periodicity in one space direction, the 
semifiow of (jlip defines a one-dimensional extended dissipative system which (nearly) satisfies 
the assumptions in Definition 12. 1[ 

Remark 2.6. The above list of examples can certainly be made longer, but all extended dissi- 
pative systems we are aware of are related somehow to a parabolic equation involving a second 
order differential operator. Higher-order systems, such as the Cahn-Hilliard equation, do not 
fit into our framework since they require a radical modification of the bound ([5]), which would 
affect our results in an essential way. 

3 Bounds on the Energy Flux 

We now begin our study of the dynamics of extended dissipative systems. Given a continuous 
semifiow {^{t))t>o on a metrizable space X satisfying the assumptions of Definition 12.11 for some 
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N G N*, we consider a trajectory u{t) = $(t)uo for which the energy density e{x,t) is uniformly 
bounded. This is always the case if our system is bounded in the sense of Definition 12.31 and as 
was observed in Remark 12.41 boundedness can often be achieved by restricting the space X to a 
suitable positively invariant subset. If we denote 

eo = sup e(x,0) < oo , /3 = sup sup6(e(x,t)) < cxd , (13) 

where e i— t- 6(e) is the nondecreasing function appearing in Definition 12.11 assumption A2 then 
implies 

\f{x,t)\^ < ^d{x,t) , xGE^, t>0. (14) 

Using only ([T]), (fT4|l . and the positivity of e{x,t) and d{x,t), we shall derive a universal bound 
on the total energy flux through a given hypersurface in R-'^ during the time interval [0,T]. 

We first consider the one-dimensional case = 1, where our hypersurface is reduced to a 
single point. Our main result in this case is : 

Proposition 3.1. Assume that N = 1, and let u{t) = ^{t)uQ be a trajectory for which the 
energy density e{x,t) satisfies p3p . Then for any x G M and any T > 0, we have 



T 



f{x,t)dt 







< VPTeo ■ (15) 



Proof. Without loss of generality, we can assume that /3 > 0. Given any T > 0, we introduce 
the integrated energy flux 

Fi{x,T) = [ f{x,t)dt, xeR, (16) 
Jo 

which is a continuous function of x G M. For any xq G M, we shall prove that Fi{xq,T) < 
{f3Teo)^^'^. Together with the corresponding lower bound Fi{xo,T) > — (/3Teo)^''^, which can be 
established in a similar way, this gives the desired conclusion. 

For any x > xq, the energy balance equation Q with 17 = (xq, x), Ti = 0, T2 = T implies 
Fi(x,r) = Fi(xo,T)+ / (e(y,r)-e(y,0))d2/+ / / d{y,t)dtdy . (17) 

Jxo ^ ' JxQ Jo 

Since e(y, T) > and e{y, 0) < eo, the first integral in the right-hand side is bounded from below 
by — eo(x — xo). On the other hand, using (I14p and the Cauchy-Schwarz inequality, we obtain 

£d{y,t)dt > ^£ f{y,tfdt > ^F,{y,Tf , 

for all y G [xo , x] . Thus (fT7|) implies 

Fi(x,r) > Fi(xo,r)-eo(x-xo) + -^ / Fi{y,Tfdy, x > xq . (18) 

Jxo 

We now compare Fi{x,T) to the solution F(x) of the differential equation 

F'(x) = -eo + -^Fixf , X > xo , (19) 

with initial data F{xo) = Fi{xq,T). If Fi{xq,T) > {/STeo)^^"^ , then F is strictly increasing 
and blows up at some finite point xi > xq- But (fT8|) then implies that Fi(x,T) > F{x) for all 
X G (xq, xi), which leads to a contradiction because we know that Fi(x, T) is uniformly bounded 
on [xo,xi]. Thus we must have Fi{xo,T) < {(3Teo)^^'^. □ 
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Remark 3.2. Albeit elementary, Proposition 13.11 already has interesting dynamical conse- 
quences. For instance, it immediately implies that an extended dissipative system on M cannot 
have any nontrivial time-periodic orbit with uniformly bounded energy, see [15]. Indeed, for such 
a periodic orbit, the last integral in the right-hand side of (I17p grows linearly in T as T — )• oo 
(if the interval [xo,x] is large enough to include a region where energy dissipation is nonzero), 
whereas the first integral is uniformly bounded by periodicity and the flux terms are 0{T^^'^) 
by (|15p . Thus (jl7p cannot hold for sufficiently large times. 

We next investigate the analog of Proposition 13.11 in the higher-dimensional case N > 2. 
Here we consider the energy flux through the boundary of the ball Br = {x G M.^ \ \x\ < R}, 
for various values of the radius R. We recall that the Euclidean measure of the sphere OBr is 
ujnR^^'^, where 



UN 



/•oo 

and r(a) = / t"~^e"*df, a>0. 
Jo 



T{N/2) ' 

Given R > and T > 0, we thus define the integrated flux 

F{R,T) = r [ f(x,t)-^dadt , (20) 

Jo J\x\=R Fl 

which represents the total energy entering the ball Bji through the boundary over the time 
interval [0, T] (the energy leaving the ball is of course counted negatively). 

Before stating our result, we introduce the higher-dimensional analog of the differential 
equation ()19p . which (after suitable normalization) becomes 

H'{r) + ^^H(r) = -1+H(rf , r>0. (21) 
r 

The following elementary result will be established in Section [71 

Lemma 3.3. For any N £ W the differential equation ()2ip has a unique positive solution 
■ (0, +oo) —7- (0, +oo). If N >2, this solution is strictly decreasing and satisfies 



/i7v(r) = 1 + ^^ + 0(32) r^+00, (22) 



and 



r log(l/r 
A^-2 



^ if N 



as r ^ . (23) 



if N >3 



r 

Moreover, any solution of (j2ip above hjy blows up in finite time, and any solution below 
cannot stay positive. Finally, Hn is given by the explicit formula 

KN{r) 

hN{r) = . , r > , (24) 

Kn_-. (r) 
2 

where Ky denotes the modified Bessel function as defined in J21 Section 9.6]. In particular 
hi{f) = 1 CLnd h^{r) = 1 + 1/r for all r > 0. 

We are now able to state the main result of this section : 
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Proposition 3.4. Assume that N > 2, and let u{t) = ^{t)uQ be a trajectory for which the 
energy density e{x,t) satisfies (llSh for some cq > and /3 > 0. Then, for any R> and any 
T > 0, the integrated energy flux ()20p satisfies 

where hj^ is given by ([2T 



Proof. Given T > and Rq > 0, we shall prove inequality ()25p for R = Rq. To do that, we 
consider the energy balance equation ([2]) in the spherical shell = {x G M'^ \ Rq < \x\ < R} 
over the time interval [0, T], where R > Rq. Using the notation (j20p we obtain 

F{R,T) =F{Ro,T)+ [ (e{x,T)-e{x,0))dx+ [ [ d{x,t)dxdt. (26) 

To estimate the right-hand side of (j26p . we proceed as in the proof of Proposition 13.11 We first 
observe that 

/ (e(x,r)-e(x,0))dx > -eo\n\ = -eo^^(i?^-0 . (27) 
Next, applying the Cauchy-Schwarz inequality to ()20p and using (jl4p . we obtain the inequality 

\F{R,T)\^ < ujnR^'^T f [ \f{x,t)\'^dadt < ^TujnR^'^ [ [ d{x,t)dadt, 

Jo J\x\=R Jo J\x\=R 



which implies that 

rT r rT t-R 



(28) 



I ( d{x,t)dxdt = I I I dix,t)dadrdt > [ ^^^^^^ dr. 

Jo Jn Jo Jro J\x\=r PTujN Jro r^^ 

Thus, combining with ([27|) and ([25]) . we obtain 

F{R,T) > F(i?o,T)-eo^^(i?^-0 + -^ ^ ^^^i^ dr , R > Ro . (29) 

Now, arguing as in the proof of Proposition 13.11 we compare F(R,T) to the solution of the 
differential equation 

F\R) = -coc.nR''-' + J^^^ ' ^ > ^0 , (30) 
with initial data F{Rq) = F{Rq,T). To eliminate all parameters in (I30p . we set 

;£gL = V^ff(fi^, R>R„, 

SO that H satisfies the normalized equation (12ip . By Lemma 13.31 if 



FjRo) 



then the solution F of ([30]) is strictly positive for R > Rq and blows up at some finite point 
i?* > Rq. But in view of ([29]), ([30]) we also have F{R) < F{R, T) for all R G {Rq, R^), which is 
impossible since F{R,T) is uniformly bounded for all R G [Rq,R^:]. Thus we must have 

F{Ro,T) 



which is the desired bound. □ 
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Remarks 3.5. 

1. The bound ()25p also holds for = 1 if we set wi = 2 and hi = 1. It then asserts that the 
total energy entering the segment [—R,R] over the time interval [0,T] is bounded from above 
by 2(/?reo)^/2. This is of course an immediate consequence of Proposition 13. 1[ 

2. It is also possible to obtain a lower bound on the energy flux leaving the ball Br over the 
time interval [0,T]. For instance, using the energy balance equation ^ with = Bji and 
[Ti,T2] = [0,T], we easily obtain the estimate F{R,T) > —LO]\fR^eo/N, which is uniform in T. 

3. If = 2, it follows from ((231) that 

2/3T 




\/ PTen h^i R, ——] ~ 7 ^ as T — )- +oo . (31) 

RMS') 

In view of Proposition 13.4^ it follows that, for any given R> 0, the integrated flux F{R,T) can 
grow at most sub-linearly (like T/logT) as T — )• +oo. This is enough to preclude the existence 
of nontrivial time-periodic solutions in two-dimensional extended dissipative systems, using the 
same argument as in Remark 13. 2^ see also [15j. In contrast, if > 3, we have 

V^ohN[R^) ~ (^-2)^ as T^+oo. (32) 

In that case F(R, T) may grow linearly in time as T — )• oo, which is compatible with the existence 
of nontrivial time-periodic orbits (see [15j for explicit examples). But such solutions must be 
spatially localized, because estimate (f25]) shows that the flux per unit area F{R,T)/{ujnR^~^) 
decreases like 1 /R as the radius R of the sphere increases to infinity. 

4. The right-hand side of (j25p always decreases to zero as /3 — )• 0. On the other hand, it is 
easy to verify that the asymptotics ()3ip and (I32p also hold for a fixed T > in the limit where 
Co — >■ 0. In particular, if A'^ > 3, Proposition 13.41 does not preclude the existence of nontrivial 
solutions emerging from initial data with zero energy density. 

5. As was mentioned in Remark I2.2[ 5. when solving a nonlinear PDE it is often convenient 
to use a function space where the energy density is not continuous, but only locally integrable. 
Although we do not want to address such technicalities in the present paper, it is perhaps 
instructive to see how Proposition 13.41 is modified if we only suppose that the initial energy is 
bounded in the uniformly local sense, namely 



sup / e{y, 0) dy < oo . 



eo 

\y-x\<l 



For simplicity, we still assume that (jl4p holds for some /? > 0. Arguing as in the proof of 
Proposition 13.41 we fix i? > i?o ^ 1 and we consider the energy balance equation in the spherical 
shell r2 = G \ Rq < \x\ < R} over the time interval [0, T]. The main difference comes from 
the estimate of the energy initially contained in fi. Using the definition of cq, we find 



/ e{x, 0) dx < c^eo^(i?^ - <) + d^eoWiv<~' 
Jn 



where CN^d^ are positive constants related to the optimal covering of a (large) ball or sphere in 
with balls of unit radius [11] . Inserting this estimate in ()27p and proceeding as before, we 
can show that F{R,T) > F{R) for R > Rq, where F is the solution of the ODE 

F'iR) = -eoCNCONR''-^ + ^rL^ ^RN-l ' ^ > , 
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with initial data F{Rq) = F{Ro,T) — e^dNOONRQ ^- This leads finally to the upper bound 



^l^^N-i ^ dNeo + VcNf3TeohN(R^I^^^ , R>1, (33) 
which replaces (j25]) . Note that the asymptotics as T — )• oo are still given by ([3T]) , (j32|) . 



4 Bounds on the Energy Dissipation 

As is clear from the balance equation a bound on the amount of energy entering the ball 
Br = {x £ M.^ I \x\ < R} over the time interval [0, T] implies an estimate of the energy dissipated 
in Bpt during the same time, provided the initial energy in Br is under control. In this section, 
we derive various dissipation estimates using the flux bounds established in Section [3l We also 
show that, for nonequilibrium solutions of extended dissipative systems on M or M?, energy 
dissipation must occur "almost everywhere" in space. 

4.1 Energy dissipation in fixed or increasing domains 

As in Section [3l we consider a trajectory u{t) = ^{t)uo of an extended dissipative system 
satisfying the uniform bounds ()13p for some cq > and /? > 0. Given R > and T > 0, we 
denote by F{R, T) the energy entering the ball Br (through the boundary BBr) over the time 
interval [0,r]. This quantity is defined by ([20]) for > 2, and if A^ = 1 we set F{R,T) = 
Fi{R,T) — Fi{—R,T), where Fi is given by ([16]). Using the energy balance equation ([2]) and 
Proposition 13.41 we easily obtain 



[ [ d{x,t)dxdt = F{R,T) + [ (e{x,0)dx-e{x,T))dx 

Jo J Br J Br ^ ' 



(34) 



where /iat is given by (p^ . Equivalently, if hiy{r) = NhN{r)/r, we find 



dix,t)dxdt < eoi^fiN^R^l^) +1^ . (36) 



<^ivi?^ Jo J Br 

We now investigate a few consequences of the general bound (j35p or ()36p . 

First, we fix > and consider the limit where T — )• +oo. Using the asymptotics ()23p for 
A > 2 and the fact that /ii = 1, we obtain the following result. 



Corollary 4.1. Under the assumptions of Proposition [3n\ or \3.4\ the following inequalities hold 
for any R > 0. 

1) IfN=l, 

2) IfN = 2, 

3) IfN>3, 



limsup^ /" / d{x,t)dxdt < 2y^^ . (37) 

T^oo VT Jo J Br 

limsup-^^ / / d{x,t)dxdt < Att/S . (38) 

limsup;^ /" [ d{x,t)dxdt < p{N -2)u}nR^^^ ■ (39) 
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In particular, if < 2, it follows from (j37p . (jSSp that the energy dissipation in any fixed ball 
converges to zero "on average" as time goes to infinity. Since we assumed that energy dissipation 
vanishes only on equilibria of the system (see assumption A3 in Definition I2.ip , these estimates 
will imply that the trajectory u{t) converges "on average" to the set of equilibria as t — )• oo, 
in a sense that will be specified in Section [5j Observe also that the bounds (j37p and (j38p are 
independent of the radius R of the ball, whereas (p8]l and (|39]l do not depend on the initial 
energy density cq. 

It is also instructive to estimate the energy dissipation in a ball whose radius depends on 
the observation time T. In view of (j36p . it is natural to take R = RoVT for some Rq > 0. We 
thus find : 



Corollary 4.2. Under the assumptions of Proposition [BTl\ or \3.^ the following inequality holds 
for any N ^'H* , any Rq > 0, and any T > : 

N f u ...... ^ 



dix,t)dxdt < eo(hN(^Ro^) +1^ , (40) 



where hN{r) = NhNir) /r . 

Observe that the volume of the space- time cylinder Br x [0,r] is lonN'^R^T, hence (gO]) 
implies that the energy dissipation rate d{x,t) is very small on average on Br x [0,T] if T ^> 1. 
This remark will be exploited in Section [6l on a particular example, to prove convergence to 
equilibria uniformly on large domains (whose size increases with time). 

Finally, in the two-dimensional case, it is also useful to consider the energy dissipation on 
a ball whose radius R has a slower growth than T^/^ as T — )• oo. Obvious possibilities are 
R = RqTi for 7 < 1/2, or = RqT^I'^ / \og{T). This gives the following estimates, which 
complement ([38|) . 



Corollary 4.3. Assume that N = 2. Under the assumptions of Proposition the following 
inequalities hold: 

1) If R{T) = RqT'^ for some Rq > and some 7 € [0, 1/2), then 

log(r) f ^ / 47r/3 



limsup — — — / / d{x,t) dx dt < 
r-5>oo T Jq Jbj^t) 

2) IfR{T) = RoT^/'^/log{T) for some Rq > 0, then 



1-27 



limsup^i^^^i^^ r [ d{x,t)dxdt < 27r/3 . 



4.2 Spatial distribution of energy dissipation 

In Section [3l we have seen a first way to exploit the energy relation ([TJ and the flux bound (fHp 
to derive useful information on the dynamics of the system. We now consider the problem from 
a somewhat broader perspective. Let u{t) = $(t)tto be a trajectory of an extended dissipative 
system in the sense of Definition 12. 1^ and suppose that the energy flux satisfies (I14p for some 
/3 > 0. This is the case if the function e 1— )• 6(e) in assumption A2 is bounded from above, or 
in more general situations if the energy density e{x, t) is uniformly bounded. Given R > and 
T > 0, we consider the integrated energy flux F{R,T) defined by (1201) . and we also denote 

E{R,T) = [ e{x,T)dx , 6E{R,T) = E{R,T) - E{R,0) . 
Jbr 
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where as usual Br = {x £ R | |x| < R}. With these notations, it follows from (j26p and (j28p 
that ^ ^ 

F{R, T) > F{Ro, T) + 5E{R, T) - 6E{Ro,T) + — ^ / dr , (41) 

for all R > Rq > 0. If we had equality in (j4ip . we could differentiate both sides with respect to 
i? and obtain, after omitting the T-dependence, the Riccati differential equation 

As a matter of fact, if F is the solution of (j42p with initial data F{Rq) = F{Rq,T) > 0, it follows 
from dH]) that F{R,T) > F{R) for i? > i^o, as long as F{R,T) + F{R) > 0. This comparison 
principle imposes strong constraints to the possible solutions of (HH), because (as we have already 
seen) the solutions of the Riccati equation may blow-up in finite time. Unfortunately, the 
conditions preventing a blow-up are not easy to specify in general, because the solutions of (P2]) 
cannot be written in explicit form. 

One way to proceed is to make simple assumptions on the source term 5E'{R) allowing 
to obtain an explicit solution of ()42p . to which the solution of (I4ip can then be compared. 
In Section [3l for instance, we assumed that 6E'{R) > —uj]\fR^~^eo for some cq > 0, and we 
obtained as a consequence the upper bound ([25]) . Here we use the same strategy to prove that, 
if < 2, the energy difference 6E{R, T) must be negative for most values of the radius i? > 0. 
Given T > 0, denote 

Jj, = |i? > E{R,T) > S(i?,0)| C (0,oo) . (43) 
The main result of this section is : 

Proposition 4.4. Assume that uq £ X is not an equilibrium, and that the trajectory u{t) = 
^{t)uQ satisfies (jl4p for some /3 > 0. Then for any T > we have 

dr < oo , (44) 

where Ijj. is the characteristic function of the set Jt defined in (j43p . 

Remark 4.5. Of course, the conclusion of Proposition 14.41 is interesting only if < 2. If A^ = 1 
then (|44p simply means that the Lebesgue measure of the set Jt C (0, oo) is finite. If A^ = 2 we 
have 

meas(J. n[l,r])^^^ T ^ dr < oo , 



which implies (roughly speaking) that meas( JtH [1, -R]) = o{R/ log(i2)) as ii — >■ oo. In both cases 
(f44l) shows that Jt is a very sparse subset of the half-line (0, +oo), so that E{R,T) < E{0,T) 
for most values of -R > 0. This considerably strengthens the results obtained (on a particular 
example) in [151 Section 2]. 

Proof of Proposition 14.41 Fix T > 0. We start from the energy balance equation (|34p . 
namely 

F{R,T) = 6E{R,T) + [ [ d{x,t)dxdt, R>0. (45) 

Jo J Br 

Since no G X is not an equilibrium, assumption A3 in Definition 12.11 implies that the last term 
in ([15|) is positive when R > Ri, for some (sufficiently large) Ri > 0. If Jt C (0,i?i], then 



14 



obviously (1^) holds. If this is not the case, we choose i?2 G Jt H (i?i,+oo) and ()15|) then 
implies that F{R2,T) > 0. Taking the limit i?o ~^ in (I4ip we also have 

1 Fir r)2 

F{R,T) > 6E{R,T) + —^ dr , R>0. (46) 

(3TujN Jo r^^ ^ 



We now define 

1 F{r,Tf 

The function : (0,oo) — )• is nondecreasing and J-{R) > for all R> R2- Moreover, using 
and the definition of Jt, we easily find 



T'{R) >lj,iR)^^ , R>0 



Thus, for all R > R2, we have 



^ dr < / — \4 dr = —— - —— < 



IR, r^-' - Jr, Hr? HR2) HR) - HR2) ' 

and (IHD follows. □ 



Remark 4.6. In the proof of Proposition 14.41 we used in fact a slightly stronger version of 
assumption A3 in Definition [2TT]: If u{t) = ^{t)uQ is a trajectory for which the energy dissipation 
satisfies d = on X [0, e] for some e > 0, then uq is an equilibrium. With the original 
formulation, the conclusion of Proposition 14.41 would hold only for sufficiently large T > 0. 
Actually, in all concrete examples we are aware of, both versions of assumption A3 are equivalent. 



5 Convergence to Equilibria 

So far we only considered a single trajectory u{t) = (^{t)uQ of an extended dissipative system, 
and under appropriate boundedness assumptions we established a few estimates on the transfer 
and the dissipation of energy. Now, in the spirit of Remark 13.21 we want to show that these 
results impose nontrivial restrictions on the dynamics of the whole system, at least if the space 
dimension is not larger than 2. In particular, we shall use the topology of the underlying space 
X to formulate convergence results, and to study the dynamics of the system in a neighborhood 
of a given point. 

For definiteness, we assume henceforth that our system is hounded and compact in the sense 
of Definition 12.31 As was already mentioned in Remark 12.41 boundedness is easy to achieve in 
the applications by restricting the space X to a suitable positively invariant subset on which 
the energy density is uniformly bounded, and compactness can then be obtained by equipping 
X with a sufficiently weak topology. If u{t) = <t{t)uQ is a trajectory of our system and if d{x^ t) 
is the corresponding energy dissipation rate, we denote for all ii > and all T > : 

D{R,T) = I AR{u{t))dt, where AR{uit)) = [ d{x,t)dx. (47) 

Jo J Br 

Here, as usual, Br = {x G | |x| < R}. By assumption Al in Definition 12.11 we know that 
D{R,T) depends continuously on the initial data uq in the topology of X. 

As a consequence of Corollary 14.11 we first estimate the time spent by any trajectory in a 
neighborhood of a nonequilibrium point. 
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Proposition 5.1. Consider a bounded extended dissipative system on M.^ with N < 2. If 
u & X is not an equilibrium point, then u has a neighborhood V in X such that any trajectory 
u{t) = ^{t)uo satisfies 



lim sup ■ 

T-s>oo T 



[ lviu{t))dt < oo , (48) 
Jo 



where ^'i(T) = vT, ^2(7") = log(r), and ly denotes the characteristic function ofVcX. 

Proof. We proceed as in [151 Section 5.1]. If u G X is not an equilibrium point, then assump- 
tion A3 in Definition 12.11 implies that the trajectory u{t) = ^{t)u satisfies D{R,Tq) > for some 
R > and some Tq > 0, where D{R,Tq) denotes the energy dissipation (jl7|l for the solution 
u{t). By continuity, there exists e > and a neighborhood V of n in X such that, for any uq G V, 
the solution u{t) = $(i)tto satisfies D{R,To) > e > 0, where D{R,Tq) is given by (l47l) . 

Now, let u{t) = ^(t)uo be any trajectory of our system. Using the notation (jI7|) . we have 
for all T > : 

AR{u{t))dt>- — / AR{uis))ds) dt 



hence 



lvKt))(^^^ \j,{u{s))dsjdt>^J^ lv{u{t))dt, 
[ lv{u{t)) dt < — [ °[ d{x,t)dxdt , r>0. 

Jo e Jo Jbh 

If we multiply both sides by 1/VT (if = 1) or log T/ \fT (if N = 2) and take the limit T — )• oo, 
we obtain (j48p using Corollarv 14.11 □ 

Remark 5.2. In [15J, the following weaker result was obtained for a particular system: If 
N < 2, any nonequilibrium point has a neighborhood V in A such that any trajectory u{t) 
satisfies 



hm i / Mu{t))dt = . 



This of course follows from ()48p . which gives a much more precise estimate of the fraction of 
time spent by the trajectory u{t) in the neighborhood V. 

Proposition 15.11 was obtained without any compactness assumption. If we now suppose that 
the space X is compact, we can use (j48p to prove that all trajectories converge in some sense to 
the set equilibria. Indeed, given any trajectory u{t) = ^{t)uQ, we can define the omega-limit set 



UJ = \ u £ X 



3 tn — ^ +00 such that n(t„) > u in X> C X . (49) 

It is known 119] that to is nonempty, compact, connected, fully invariant under the semiflow 
$(t), and that distx {u(t),uj) — )■ as t — )■ -|-oo. However, our assumptions do not imply that uj is 
contained in the set of equilibria. Counter-examples can indeed be constructed even for relatively 
simple systems such as the Allen-Cahn equation in one space dimension, see Example 15. 71 below. 
Motivated by the conclusion of Proposition 15. H we propose the following alternative definition : 

Definition 5.3. If u{t) = ^{t)uQ is a trajectory of a bounded and compact extended dissipative 
system on with N < 2, we define 

O = iu £ X limsup — ^ — ^ / l-\;{u{t))dt = oo for all neighborhoods V of u\ , (50) 



where ^'i(r) = VT and ^'2(r) = log(r). 
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In other words, u) is the set of points in all neighborhoods of which the trajectory u{t) spends 
a "substantial fraction of the total time". What is exactly meant by "substantial" depends on 
the space dimension A^, and is specified by the function ^i\}{T). It is clear from the definition 
that (D C w, and Proposition 15.11 implies that uj is contained in the set of equilibria of our system. 
More properties of a) are collected in our final result : 

Proposition 5.4. Let u{t) = ^{t)uQ be a trajectory of a bounded and compact extended dissi- 
pative system on with N < 2. Then the set OCX defined by (150p is nonempty, compact, 
and contained in the set of equilibria. Moreover, ifV is any neighborhood of ui in X, then 

lim sup J; ' / lv4u{t))dt < oo . (51) 



T 



Proof. We proceed as in [T5l, Section 5.2]. We first observe that, if T C X is compact and does 
not intersect Co, then there exists a neighborhood V of T such that 

lim sup ]} ' \ lr('u(0)dt < lim sup J; ' \ \v{u{t))dt < oo . (52) 

Indeed, this property holds by definition if F = {ui} for some ui ^ Co, and the general case follows 
by a finite covering argument. Now, if we take for T the closure of the trajectory {u{t) \t > 0}, 
then r is compact and T^^ Jq lr{u{t)) dt = 1 for all T > 0, which is incompatible with ()52p . 
Thus we must have T Co ^ 0, hence in particular w 7^ 0. Moreover, it is clear from the 
definition that Co is closed in X and contained in T, hence Co is compact. On the other hand. 
Proposition 15.11 precisely means that Co is contained in the set of equilibria. Finally, if V is any 
open neighborhood of Co in X, then F n is compact and does not intersect Co, hence by ()52p 

^'Ar(T) /■'^ ^'Ar(T) f'^ 

limsup — — — / lv{u{t))dt = limsup — — — / l\>cf-^j-{u{t)) dt < 00 , 
r-5.oo T Jq T-^OO T Jq 

which proves ()5ip . □ 

Remark 5.5. Since Co consists of equilibria, it is obvious that ^(t)Co = Co for all t > 0. In 
fact, for any relatively compact trajectory of a continuous semifiow on a metrizable space X, 
one can prove that the set Co defined by ()50p is nonempty, compact, and fully invariant, see \15\ 
Proposition 5.4]. These properties are therefore independent of the gradient structure. On the 
other hand, the set Co (unlike lo) is not connected in general, as can be seen from Example 15.71 

Remark 5.6. Instead of Co, the following set was defined in [15] (for a particular system) : 
00 = iu € X limsup — / l\i{u{t)) dt > for all neighborhoods V of "ul . 

^ T-s>oo T Jq ) 



Clearly Co <Z uj, hence Proposition 15.41 implies that Co is contained in the set of equilibria, as was 
proved in jl5l Proposition 5.4]. It is also known that 7^ 0, which implies that w 7^ 0. 

Example 5.7. Under the assumptions of Proposition 15.41 it is not a priori obvious that the 
usual omega-limit set (j49p is not necessarily contained in the set of equilibria. In this respect, the 
following example is instructive. We consider the one-dimensional reaction-diffusion equation 
^ with the double- well potential V{u) = \{l - u^ f : 

dtu = dlu + u-u^ , X G M , t > . (53) 
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This system has three constant steady states : uq = (which is unstable), and u± = ±1 (which 
are stable). In addition, there is the "kink" solution 

ipix) = tanh(x/\/2) , x G M , (54) 

which connects U- at x = — oo to at x = +oo. It can be shown that u± and the translates 
of ±'0 are the only stable steady states of ([53]) in the space of bounded solutions. 

Interesting nonequilibrium solutions of ([53|) can be constructed by gluing widely separated 
kinks. For instance, if a ^ 1, the function 

Va{x) = ip{x - a) - '4){x + a) + 1 , x G R , 

describes the superposition of a kink located near x = a and an "anti-kink" —tjj near x = —a. 
This is not an equilibrium of ()53p . but it can be shown that the solution of ()53p with initial 
data Va stays very close to ya{t) for later times, provided the parameter a evolves according to 
the exponential law a ~ — ci exp(— C2a), for some ci,C2 > 0, see e.g. [8]. This approximation 
property remains valid as long as both kinks are widely separated, but when they get close to 
each other they "annihilate" and the solution converges uniformly to 1 as t — ?■ +oo. 

Using these results and a general procedure that can be found e.g. in [13], one can show 
that there exists a unique eternal solution : M x M — )■ M of (j53p such that ti^(0, 0) =0 and 

sup|nw,(x,t) — Va(t){x)\ > , where a{t) — )• +oo as i — )• — oo . (55) 

In fact, one has a{t) ~ C2 ""^ log(|t|) as t — t- —00. This solution converges uniformly to n_|_ = 1 
as t — 7- +00, and uniformly on compact sets to ti„ = — 1 as t — >• —00. If 7Ioc(R) denotes the 
topology of uniform convergence on compact sets of M, it follows that u^{t) realizes a heteroclinic 
connection from n_ to ti+ through the symmetric collapse of a pair of kinks coming from infinity. 

Now, using an idea taken from |10j . we consider the solution u : M x M+ — )■ M of ()53p with 
initial data uq satisfying 

no(x) = (-1)"+! if bn < \x\ < bn+l , (56) 

where (6n)neN is a strictly increasing sequence satisfying 60 = and bn+i 3> bn for all n G N. 
Under the evolution of the parabolic equation (|53p . the discontinuities of the initial data are 
rapidly smeared out, and replaced by smooth interfaces of the form (I54p . the positions of which 
slowly move according to exponential law specified above. In particular, since bi ^ bo, the pair 
of kinks closest to the origin will hardly feel the presence of the other kinks, and will therefore 
evolve in time like the solution of ()53p with initial data V^q. Once the first pair has disappeared, 
we are essentially back to the original configuration, with a central pair of kinks that is now 
located near ±61. This pair evolves on a much slower time scale, but will eventually come close 
to the origin and annihilate, and the same process will continue forever since we started with 
infinitely many kinks. Such a coarsening dynamics was studied for instance in |121 125j. 
These heuristic considerations lead to the following reasonable conjecture : 

Conjecture 5.8. Let n : M x M_(_ — )• M 6e the solution of (j53p with initial data (|56p . Then the 
omega-limit set of the trajectory {u{t))t>o in the topology Tiod^) is 

oj = {u+ , n_} U {u^{t) 1 1 G M} U {-u^{t) 1 1 G M} , 

where u± = ±1 and is the eternal solution of (j53|) defined by ([55]) . 
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If this conjecture is true, then u) consists of two equihbria u± and two heterochnc connections 
between them. Thus w is a heterochnic loop, which is not entirely contained in the set of equilib- 
ria. In contrast, for the same solution, the modified omega-limit set introduced in Definition 15.31 
satisfies u = {«+ ,"«_}, hence is contained in the set of equilibria. Note that Proposition 15.11 
implies that the number of annihilations of pairs of kinks that can occur in the time interval 
[0, T] is bounded by C^T for large T. 



6 The Vorticity Equation in an Infinite Cylinder 

In this section we analyze in some detail an interesting example which does not fit exactly into 
the framework of Definition 12.11 but can nevertheless be studied using the techniques developed 
in Sections [3] to El We consider the incompressible Navier-Stokes equation in the infinite cylinder 
O = M X T, where T = M/Z. Points in Vt are denoted by x = (xi,j;2)) where xi S M is the 
"horizontal" and X2 G T the "vertical" variable. Our system reads 

dtu + {u ■ V)u = An — Vp , divu = , (57) 

where u : fix R_|_ denotes the velocity field and p : fix M_|_ — )• M the pressure field. For each 

t > 0, both quantities u{xi, X2,t), p{xi,X2,t) are assumed to be bounded in Q and 1-periodic 
with respect to X2- Since u is divergence free, we have in particular 

di / Ui{xi,X2,t)dx2 = / diUi{xi,X2,t)dx2 = - d2U2{xi,X2,t)dx2 = , 

Jt Jt Jt 

hence the vertical average of the horizontal speed, which we denote by (ni), does not depend 
on the horizontal variable xi. As is explained for instance in [2], it then follows from ()57|) that 
dt{ui) = 0, so that {ui) is a constant which can be set to zero by an appropriate Galilean 
transformation. Thus, without loss of generality, we can assume that 

n(xi,X2,i) = f .V.') + f^'!^''''''!ll ' {xi.X2)(^fl, t>0, (58) 

\m{xi,t)J \U2{xi,X2,t)J 

where m = (^2). By construction, we then have {ui) = {U2) = 0. 

In addition to (jSTj) . we shall study the evolution equation for the vorticity uj = diU2 — 92^1- 
In view of (I58p . we have 

Uj{xi,X2,t) = dim{xi,t) -\- U){xi,X2,t) , {xi,X2)£fl, t>0, 

where dim = (u) and co = diU2 — d2Ui- It is important to realize here that, since we want to 
consider solutions of ([57]) which do not necessarily decay to zero as — )• 00, the velocity field 
u is not entirely determined by the vorticity cj. More precisely, one can show that the oscillating 
part u of the velocity field is given by a Biot-Savart formula : 

u{xi,X2,t) = / V-^K{xi-yi,X2-y2)ioiyi,y2,t)dy2dyi , (59) 
where V"*- = {—82,81) and 

K{xi,X2) = — log(^2cosh(27rxi) - 2 cos(27r2;2) j — , ixi,X2) G fl , (60) 

see [2j. However, the vertical average m = {U2) cannot be completely expressed in terms of the 
vorticity, and we only know that dim = (cj). The following estimates will be useful. 
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Lemma 6.1. There exists a constant Ci > such that, for any u G L°^{Q), the velocity field u 
defined by (159]) satisfies 

W'^Wi'^^in) < Ci||^llL°°(n) • (61) 
Moreover, we have ui = —d2V for some v £ L°°{Q), and there exists C2 > such that 



lbllL°°(n) < C2esssup(/' ^2)^ dx2^ < C2\\uj\\L^(n) • (62) 

Proof. Here and below we denote by * the convolution on 17 (considered as an additive group). 
As is easily verified, the function K defined by (1601) satisfies G L}-(Vl), hence 

II^*IIl°°(0) = l|V"^i^ * tl'llLoo(Q) < ||Vi^||/^i(f^)||d)||ioo(f^) < 2||Vi^||ii(f^)||a;||ioo(Q) . 

Similarly, one can check that K G L^(0) and 

Mq := / sup \K{xi,X2)\dxi < 00 . 
Jr x2eT 

Thus, if we define v = K *uj, we have ui = —d2V by (j59]) and a standard calculation shows that 
llt'llioo/m < Mo||K||2,i(Q) esssup / cD(xi,X2)^dX2 , 

which is the desired result. □ 

Instead of the Navier-Stokes equation ()57p . we now consider the evolution system satisfied 
by the vorticity w and the average speed m = (u2)- As in [2j we obtain 

I dtm + di{uiU2) = dfm, xi G M , ^^^^ 
I dtuj + u-Vuj = AU! , {xi,X2)£^. 



Here it is understood that u is given by ([58]) . where ui,U2 are obtained from u via (l59]) . Note 
that system (j63p is somewhat redundant, because the horizontal derivative of the first equation 
is the vertical average of the second one, but as is explained above it is not possible to get rid 
completely of the first equation. Given a solution of ()63p . we define for all xi G M and t > 0: 

e{xi,t) = I: I uj{xi,X2,t)'^ dx2 , 

f{xi,t) = 1^1 {diu? - uiu?^{xi,X2,t)dx2 , (64) 

d{xi,t) = / \Vuj{xi,X2,t)\'^ dx2 ■ 
Jt 

In agreement with the general terminology used in this paper, we shall call e the energy density, 
/ the energy flux, and d the energy dissipation rate, although the term "enstrophy" would 
certainly be more appropriate than "energy" in the present context. Using (I63p . it is easy to 
verify that the quantities (j64p satisfy dte = dif — d, which is the one-dimensional version of ([T]). 
On the other hand, if d = 0, then certainly dtUJ = and Co = u — (w) = 0. Then n = by ([59]) . 
and since dfm = di (uj) = it follows from (j63p that dtm = too. Thus d = only for equilibria 
of system (I63p . Finally, we have the following estimate for the energy flux : 
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Lemma 6.2. There exists a constant C3 > such that 

< Csfl + sup e(yi) )e(xi)(i(xi) , for all xi G 



(65) 



Proof. We fix xi G M and consider both terms in (|64p separately. First, using the Cauchy- 
Schwarz inequality, we easily find 



dioj"^ dx5 



T 



ujdico dx5 



T 



< (2e)V2di/2 



On the other hand, since ui = ui = —d2V by Lemma l6.ll we have 

1 
2 



UlUJ dX2 



T 



{d2v)uj'^ dX2 



T 



VUjd2UJ dX2 



T 



hence using ([62|) we conclude 



- / UlU) dX2 



< V 



Uoo(^)(2e)V2di/2 < 2C2Supe(yi)V2eV2di/2 



Combining both estimates we obtain (j65p . 

The Cauchy problem for Eq. (j57p is globally well-posed in the Banach space 



□ 



X = {tz G CU^)' 



divu = 



equipped with the L°° norm, see [21 [T71 [Ml 127] . If u{t) is the solution of (f57|) with initial data 
uo G X, it is known that ||n(t)||Loo cannot grow faster than t^^'^ as t — )• 00 (see [2] and ([66]) below), 
but otherwise we have no information on the long-time behavior of the solution. In particular, 
uniform boundedness is an open problem, which we hope to address in a future work. Here our 
goal is to obtain some information on the associated vorticity u{t). Without loss of generality, 
we assume that the initial vorticity ujq = curluo is bounded, and we denote M = \\ujo\\l°°- Since 
uj{t) evolves according to the advection-diffusion equation (j63p . the maximum principle implies 
that ||a;(t)||Loo < M for ah t > 0. It then follows from that \\u{t)\\Loo < CiM for ah t > 0, 
so that the oscillating part of the velocity is under control. On the other hand, if we apply 
Duhamel's formula to the first equation in (|63p . we obtain 



m{t) 



e ^ mo 



die'~*-'^^-ui{s)u2{s)ds 



t > 



where tuq = m(0) is the vertical average of the vertical initial speed uq. The uniform bound on 
u{t) thus implies 



mi 



< llm-olli" 



\u[s] 



ds < ||mo||L°° + 



t > 



vr 



(66) 



/o i/7r(t - s_ 

hence ||ii(t)||Loo grows at most like t^/^, as already announced. 

Under our assumptions, the energy density and the energy flux defined by ([6l[) satisfy the 
following uniform bounds e{xi,t) < cq and /(xi,t)^ < (5d{xi,t), where 



eo 



1|| l|2 

-||a;o||L°° 



and 



/? = C73eo(l + eo) , 



(67) 



see Lemma 16.21 Thus we are exactly in position to apply the results of Sections [3] and [H In 
particular, using ([35]) with iV = 1, we obtain : 
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Proposition 6.3. // the initial vorticity is bounded, the solution of (j63p satisfies, for allT > 
and all R > 0, 

[ [ [ |Va;(xi,X2,t)pd,T2d.Tidt < 2V/3Teo + 2Reo , 
Jo J-rJj 



where eo,/3 are given by ([6 

Proceeding as in Section [5l one can then use ([68]) to show that the vorticity uj{xi,X2,t) 
converges uniformly on compact subdomains toward the set of equihbria 



£ = ^^weC^^{n) Vt« = 0, |w| < IIwoIIl-} . 



More precisely, adapting Proposition 15.41 to the present situation, we infer that, if V is any 
neighborhood of C X in the topology of C^(r2), the fraction of the time interval [0,T] spent 
by the trajectory uj{t) outside V does not grow faster than CT^^"^ as T — )• oo. This already gives 
valuable information on the solutions of (j63p , but combining (j68p with further a priori estimates 
one can obtain a stronger and more precise conclusion. In what follows, we assume that the 
solution of (f63]l under consideration satisfies 

sup / \d'l2(^{xi,X2,t)\'^ dx2 < Ml , and sup |m(xi,t)| < M2(l + t)^ , (69) 

for ah t > 0, where Mi, M2 > and /3 G [0,1/2]. The first estimate in ^ is verifed for all 
i > 1 by any bounded solution of the vorticity equation, due to parabolic regularization, and 
the second estimate with ^ = 1/2 is just ([66]) . In fact, it is possible to show that ([69]) always 
holds with j3 = 1/6, see [16], but it is still unclear whether all solutions of (j63|) with bounded 
initial data satisfy (I69p with /3 = 0. 

Proposition 6.4. Consider a solution of ()63p with bounded initial data such that ()69p holds 
for some j3 E [0, 1/2]. If f3 < a < 1/2, there exists a constant Kq > such that, for any T > 1, 

measjtG [0,T] sup sup |w(xi, X2, t)| > 7R^^] < KoT^^+^Z^ . (70) 

Remark 6.5. Estimate (|70p is useful especially when f3 < a < 1/2. It then shows that the 
vorticity uj{xi,X2,t) converges uniformly to zero on subdomains of size 0(t(""'"^^)/^), at a rate 
comparable to except for possible excursions whose probability density decays roughly 

like as i qc The fact that w converges to zero, and not to a nonzero constant w, can 

be understood as follows. If u{xi,X2,t) = t(; / on a sufficiently large spatial domain, then 
dim = (uj) = w for xi in a large interval, and this is compatible with the assumed upper bound 
(f69l) only if w is small enough. Thus, in the particular case of equation (f63l) with bounded initial 
data for uj and m, the omega-limit set "on average" as defined in Remark 1 5 . 6 1 consist s of a single 
point. 

Proof of Proposition 16.41 Applying (|68l) with R = VT, we see that there exists C4 > such 
that 

/ / / |Vw(xi,X2,t)|^dx2dxidt < dVr , (71) 
Jo J-VtJt 



I-Vt 

for all T > 1. Given a G [0, 1/2], we define 



JaiT) = {te[o,T] 



J ^J^N^{xi,X2,t)\'^dx2dxi > ^1 C [0, 



T] . 
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It follows from ([71]) that meas{Ja{T)) < C4r°+^/2, for all T > 1. Our goal is to give a uniform 
bound on the vorticity uj{xi,X2,t) on a large spatial domain for all t G [0,T] \ Ja{T). 

We observe that \uj{xi, X2,t)\ < \g{xi,t)\ + h{xi,t) for all (xi,X2) G and all t G [0,T], 
where 

g{xi,t) = f uj{xi,X2,t)dx2 , and h{xi,t) = ( f \d2Uj{xi,X2,t)\'^ dx2) 
We first bound the average g = {uj). If L < \/T and t G [0, T] \ Ja{T), we have 

[ \dig{xi,t)\^ dxi < [ [ \dMxi,X2,t)\^dx2dxi < ^ . (72) 

J~L J-L JT ^ 

Furthermore, we know that g = dim, where m{xi,t) satisfies (j69p for some /3 G [0,1/2]. Using 
(j69p . (|72p and Lemma 16.61 below, we thus obtain 
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sup |5(xi,t)| < ^ + ^-^ , tG [o,r]\ j„(r) , 

for some C5 > 0. If we now choose L = 7"("+2/3)/3 < T"^/^^ -^g arrive at 

sup{|5(xi,t)|||xi|<r("+2/3)/3| < tG[o,r]\j„(r), (73) 

for some Ce > 0. 

On the other hand, we know that J_^h{xi,tf dxi < T"" when t G [0,T] \ Ja(T). In 
addition, it follows from (j69|) that 



|5i/i(xi,t)| < (y |5^2'^(^i'2;2,t)pdx2y''^ < Ml , xiGM, t G [0, T] . 
Thus Lemma 16.71 below implies that 

sup |/i(xi,t)| < CmaxI ,„ , , ,r, /r, ) = ;o ; 

tG [o,r]\^(r) , (74) 



for some Cy > 0. Combining ([73]) . (j7i|) we obtain 

Cg C7 Cg 

sup sup |u;(xi,X2,t)| < sup b(xi,t)|+ sup |/i(xi,t)| < J + ^ < (^^73 ' 

for all t G [0, T] \ J„(r). Since L = T^<^+^^)/'^ and meas(J„(r)) < (:74T"+^/2^ this gives □ 

Finally, we state and prove two elementary interpolation lemmas which were used in the 
argument above. 

Lemma 6.6. Assume that g G C^([0,L]) satisfies 

I g'{x)'^dx < e , and / 5f(x) dx 
Jo Jo 

Then sup \g{x)\ < ^ + (Le)^/^ 



L 

< M 
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Proof. We decompose g{x) = g + h{x), where g = L ^ g{y) dy. Since h has zero mean over 
[0, L], there exists xq £ [0, L] such that h{xo) = 0. For all x G [0, L], we thus have 



\h{x)\= / h'{y)dy <\x-xo\'/^( h'{yfdy) < {Lef^ . 

Since |^| < M/L by assumption, we obtain the desired result. □ 
Lemma 6.7. Assume that h G C^{[Q,L]) satisfies 

I h{x)'^ dx < e , and sup |^'(2;)| < M . 

Jo 0<x<L 
Then sup \h{x)\ < max 



0<x<L 



Proof. If xq € [0,L] is a point where is maximal, we have 

> ||/i||L°° -M|x-xo| , xG[0,L], 
where ||^||l°° = sup{|/i(x)| | < x < L}. By straightforward calculations, we thus find 



M 



This gives the desired result. □ 

7 Appendix : Proof of Lemma 13.31 

In this final section, we study the positive solutions of the ordinary differential equation 

h'{r) + ^^h{r) = h{rf - 1 , r > , (75) 
r 

and we prove Lemma 13. 3i All arguments are quite standard, and are reproduced here for the 
reader's convenience. Although the unique positive solution of (j75p is given by an explicit 
formula which can be found using a Cole-Hopf transformation, we find it more instructive to 
prove the first part of Lemma 13.31 including the asymptotics (j23p and ()24p , without using this 
explicit representation, which will be derived only at the end. We proceed in several steps : 

1. Construction of the stable manifold. The nonautonomous ODE (j75p has an asymptotic 
equilibrium /i = 1 at r = +oo, with a one-dimensional stable manifold which contains precisely 
the solution /ijv we are looking for. To construct the stable manifold, we set 



and obtain for g the ODE 



h{r) = 1 + ^ + g{r) , (76) 
2r 



/M = 2»M+»M=-<^^g^. (77) 



As is easily seen, any solution of ([771) that stays bounded as r — )• +oo satisfies the integral 
equation 



C30 
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Now, fix R e (0,1) and take ro > large enough so that (iV-l)(7V-3) < 4Rr^. It is then 
straightforward to verify that the right-hand side of ()78p defines a strict contraction in the 
closed ball 

Br,{R) = {ge C70([ro, +00)) sup \g{r)\ < r} , 

hence has a unique fixed point G Br^iR) which, by construction, is a solution of ([77]) for 
r > rg. Since R can be taken arbitrarily small (at the expense of choosing rg large enough), it 
is clear that 5Ar(r) — >• as r — >• 00. Thus defining 

N-l 

hN{r) = 1 + ^— + gN{r) , r > ro , (79) 
2r 

we see that satisfies ([75]) and hiy{r) — )• 1 as r — >• cx). By construction /jat is the unique 
solution of ([75]) that stays bounded as r — )• +00. 

Remark 7.1. If = 1 or = 3, it is clear from ([78]) that ^'at = 0, so that hi(r) = 1 and 
hs{r) = 1 + 1/r. 



2. Asymptotic behavior as r ^ +00. Using (I78p . we easily find 

g^(r) = i^^Z^^^Z^ + , as r ^ +00 . 

Thus ([^^ holds, and in view of ([77]) we also have g'j^{r) = 0{r^'^) as r — )• +00, so that 

If iV > 2, this shows that h!^{r) < for r > sufficiently large. 

3. Global monotonicity. We assume from now on that N > 2. Solving (|75p backwards, we 
construct (for some r^, > 0) a maximal solution : (r=i,,-|-oo) — t- M which coincides with ([79p 
for r > tq. We claim that h!^{r) < for all r > r^,. Indeed, assume on the contrary that there 
exists f > r^: such that /i^(f) = and h'^(r) < for all r > f. Then 

h%{f) = (2/i7v(r) - ^^)/i'iv(^) + ^^/i7v(f) > , 

because h'j^{f) = and hN{f) > 1. This implies that h'j^{r) > for r > f close enough to f, in 
contradiction with the definition of f. Thus h'j^{r) > for all r > r^, and using ([75]) we deduce 



N-l / (7X^-1)2 
K /i7v(r) < ^ + ^1 + ^^^ forall r>r, . (80) 

This estimate shows in particular that hjy cannot blow up at a finite point r > 0, hence we 
necessarily have r^, = and h'^{r) < for all r > 0. 

4. Asymptotic behavior as r — >• 0. Setting /a? = l//iAr we obtain the ODE 

fNi^) = ^Mr) + fN{rf - 1 , r > , 
which is very similar to ([75p . In particular, we have 

(-) (l-/7v(s)2)ds, 0<r<l. (81) 
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A direct study of ([8T]) gives the following asymptotic expansion as r — )• 



f2{r) = rlog — \- Cr + o(^r'^ (log ^ for some C G 



whereas /3(r) = r/(l + r) and 



/4(r) = ^+o(r3logi) , /^(r) = ^ + 0(r3) if AT > 5 . 



Since = 1/ /n, this proves (p3|) . 

5. Uniqueness and threshold behavior. We first study the solutions of ()75|) that lie above h^- 
Assume that h is a solution of (j75p such that h{ri) > h]\f{ri) for some ri > 0. In particular, we 
have that h{r) > h]\j[r) > 1 for all r > ri. If /i'(r) < for all r > ri, then /i(r) converges to some 
limit ^ > 1 as r — )• +oo, and using (|75p we easily see that £^ = 1, hence 1 = 1. This implies that 
the function g{r) defined by (I76p is small for sufficiently large r > and satisfies the integral 
equation ([78]) . hence coincides with 5'Ar(r), which is of course impossible since g{ri) > gN{fi)- 
Thus there must exist r2 > ri such that h{r2) > h]\}{r2) and h'(r2) > 0. If we now choose 
r3 > r2 so that h'(r) > for all r G [r2,?'3], we have on that interval h{r) > (A— l)/(2r), hence 

/ A-l\ A-1 
h"(r) = h'ir) 2h(r) H 7^h(r) > , r2 < r < r^ . 

This argument shows that h{r) is convex for r > r2, and blows up at some finite point r* > r2. 
Indeed, if h{r) was defined for all r > r2, the convexity would imply that h{r) — t- +oo as 
r — )• +c«, so that h would satisfy the differential inequality h'[r) > ^/i(r)^ for all sufficiently 
large r, which is impossible because this inequality has no global positive solutions. 

We next consider solutions of ([75]) that lie below hj\f. Assume that /i is a solution of ([75|) such 
that < h{ri) < /lAf(ri) for some ri > 0, so that h{r) < h]\}{r) for all r > ri. If h[r) > for all 
r > ri, we have h'{r) < for all r > ri, hence h{r) converges to some limit i G [0, 1] as r — )■ +oo. 
But the same arguments as above show that i'^ = 1 and £ ^ 1, which is a contradiction. So the 
solution h{r) must necessarily change sign for r > ri. 

It follows in particular that h]\f is the unique positive solution of ([75]) that is defined for all 
r > 0. 

6. Explicit representation. Let u{r) = exp(— hj\[{s) ds). Then u solves the second order linear 
ODE 

// A - 1 , 

u"(r) H u'(r) = uir) , r > , 

r 

and u{r) decays exponentially to zero as r — t- +oo. Setting u = 1 — A/2 and n(r) = r"v{r), we 
obtain for v the differential equation 

r'^v" {r) + rf'(r) — (r^ + z^^)?j(r) = , r > , 

which defines the modified Bessel functions, see IJ^ Eq. 9.6.1]. Since v{r) decays exponentially 
as r — )• +00, we must have v{r) = CKy{r) for some C > 0, see [P, Section 9.7]. Thus u{r) = 
Cr^Kj^ir), and using P Eq. 9.6.28] we also find u'{r) = -Cr"" Ky^i{r). Since K^y{r) = Ky{r) 
by [U Eq. 9.6.6], we conclude that 

u(r) Ky{r) Aiv_,(r) 

2 
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which proves ()24p . The proof of Lemma 13.31 is now complete. 



□ 
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